Effects of atomic short-range order on the properties of perovskite alloys in their 

morphotropic phase boundary 
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The effects of atomic short-range order on the properties of Pb(Zri_ x Ti I )03 alloy in its mor- 
photropic phase boundary (MPB) are predicted by combining first-principles-based methods and 
annealing techniques. Clustering is found to lead to a compositional expansion of this boundary, 
while the association of unlike atoms yields a contraction of this region. Atomic short-range order 
can thus drastically affect properties of perovskite alloys in their MPB, by inducing phase transitions. 
Microscopic mechanisms responsible for these effects are revealed and discussed. 

PACS numbers: 77.84.Dy,81.30.Bx,77.80.Bh,77.22.Ej 
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Complex insulating perovskite alloys, with the general 
formula (A'A"..)(B'B"...)03, are of great interest for a 
variety of device applications because of their anoma- 
lously large electromechanical responses 0, 0, . Exam- 
ples of such applications include piezoelectric transduc- 
ers, actuators, as well as, dielectrics for microelectronics 
and wireless communication. 

Interestingly, many perovskite solid so- 
lutions, e.g., Pb(Zr 1 _ :E Ti K )03 (PZT), 
[Pb(Zn 1/3 Nb 2/ 3)0 3 ]i_ a; [PbTi03] :i; (PZN-PT), 
[Pb(Mg 1/3 Nb 2/ 3)0 3 ]i_ a; [PbTi03] a; (PMN-PT) and 
[Pb(Sc 1/2 Nb 1/2 )0 3 ]i_ x [PbTi03] :c (PSN-PT), exhibit 
their largest electromechanical responses for compo- 
sitions lying within the so-called morphotropic phase 
boundary (MPB). This area was thought for more 
than 50 years to dis ■continuously separate, in the 
temperature-composition plane, a rhombohedral (R) 
ferroelectric phase exhibiting an electric polarization 
along a (111) direction from a tetragonal (T) ferroelec- 
tric structure having a polarization pointing along a 
(001) direction. The recent discovery of a monoclinic 
ferroelectric phase in the MPB of PZT has completely 
changed this long-time accepted picture 0, since this 
new phase acts as a structural bridge between the R and 
T phases 0, Q. Furthermore, the polarization in this 
low-symmetry phase continuously rotates when varying 
the composition 0, 0| , explaining why the MPB is the 
region of choice for optimization of piezoelectric and 
dielectric responses H, IE 0> El These recent findings 
have led to a flurry of investigations aimed at better 
understanding and characterizing the properties of 
the MPB in various perovskite alloys. In particular, 
other low-symmetry phases have also been subsequently 
discovered in PZN-PTJPMN-PT and PSN-PT (see, for 
instance, Refs. H [H E H Q 03 ) . 

One remaining mystery of the MPB is about its com- 
positional width. More precisely, rather different widths 
have been reported for the same solid solution, depend- 
ing on the growth conditions Chemical short-range 
ordering (SRO) is often invoked for this variance between 



different measurements [l^ | , since no long-range ordering 
occurs (to our knowledge) in the mixed sublattice of PZT 
and PMN-PT near their MPB. However, why and how 
SRO would affect the morphotropic phase boundary are 
two unresolved questions. One possible reason for this 
lack of knowledge is that the characterization of SRO 
is challenging and can only be accomplished via non- 
conventional experiment [17]]. Another plausible reason 
is that mimicking these effects — via the use of compu- 
tational schemes — requires high accuracy and handling 
of large supercells, which are two conditions that are not 
simultaneously met by cither usual first-principles tech- 
niques or semiempirical approaches. 

In this Letter, we take advantage of the accuracy, the 
possibility of using large supercells, and the microscopic 
insight provided by the first-principles-based approach 
of Ref. 5] to study the effect of SRO on the physical 
properties of PZT in its MPB. The use of this technique 

(1) proves that SRO does have a drastic effect on these 
properties, especially in the Ti-poor region of the MPB; 

(2) further indicates that short-range association of like- 
atoms has an opposite effect than short-range association 
of unlike- atoms on the compositions delimiting the MPB; 

(3) reveals the microscopic mechanisms responsible for 
these features. 

Short-range ordering in any A(B' 1 _ 2 .B")03 perovskite 
alloy can be characterized by the so-called Cowley pa- 
rameters, defined as [l8|: 



aj(x) = 1 



P(j) 



(1) 



where P(j) is the probability of finding a B' (respec- 
tively B") atom being the jth nearest neighbor of a B" 
(respectively B') atom in the mixed B-sublattice. Then, 
otj(x) = for all j's represents a truly disordered alloy 
while otj(x) > is associated with clustering , that is jth 
nearest neighbors are preferentially of the same atomic 
kind. Conversely, a configuration for which ctj(x) < 
corresponds to anticlustering , i.e. to a situation in which 
jth nearest neighbors are preferentially of a different 
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atomic kind. Note that an increase in the magnitude 
of ctj is indicative of an increase in the strength of SRO. 
We focus here on relatively small ordering deviations oc- 
curring over the first nearest neighbors shell from the 
perfectly disordered Pb(Zri_ a; Ti a ;)03 alloys. Practically 
we limit ourselves to the range —0.3 < otx < 0.3, while 
always stipulating that a.i = 0*3 = 0. We use a simu- 
lated annealing technique that allows us to deter- 
mine the 12 x 12 x 12 atomic configuration {a{\ — where 
(7;=+l or —1 corresponds to the presence of a Zr or Ti 
atom, respectively at the B-sublattice site % of PZT - 
that best satisfies the pre-fixed Cowley parameters. Once 
the desired atomic configuration is generated, it is used 
for input in Monte-Carlo simulations utilizing the first- 
principles-derived alloy effective Hamiltonian proposed in 
Refs. • All the simulations conducted in the present 
study are done at 50 K. The outputs of the Monte-Carlo 
procedure are the set of polar local soft modes {ujj, 
where i indexes the different 5-atom perovskite cells |20j , 
and the homogeneous strain tensor. The (u) supercell 
average of the {ui} modes is directly proportional to the 
macroscopic electrical polarization. The {u^} supercell 
distribution is also used to compute the (k) coefficients 
defined as: 

/ M (k) = y 8^exp(ik-R i K M , (2) 

i 

where k is a vector in the first Brillouin zone of the sim- 
ple perovskite structure and Ri is the lattice vector as- 
sociated with cell i. /1 identifies the cartesian coordinate 
(the x, y, and z axes are chosen along the pseudo-cubic 
[100], [010], and [001] directions, respectively) of the lo- 
cal modes u^. (3 is a normalization coefficient yielding 
J2k l/^(k)| 2 = 1- It i s important to realize that a value 
close to 1 for |/ M (k)| 2 at k = — i.e., for the T-point 
— corresponds to the homogeneous situation where the 
/i-component of is nearly independent of the cell i. 
On the other hand, the inhomogeneous case in which the 
local polarization fluctuates from its averaged value be- 
tween different cells yields a [/^(k)! 2 that is smaller than 
1 for k = 0. The larger this real-space fluctuation is, the 
smaller |/ M (k)| 2 at T. 

Figure 1 shows the composition dependence of the x, y, 
and z Cartesian coordinates ((u x ), (u y ), and (u z )) of (u) 
as a function of composition in disordered materials. As 
already indicated in previous studies 0, , our effective 
Hamiltonian approach successfully predicts the existence 
of three ferroelectric phases for random Pb(Zri_ a; Ti a ;)03 
solid solutions near its MPB: a tetragonal (T) phase for 
x >48.7%, a rhombohedral (R) phase for x <47.7%, and 
the recently discovered monoclinic - and so-called Ma 
[24I - phase for the composition range Aa;=1.0% in be- 
tween. The polarization is parallel to the pseudo-cubic 
[001], [111], or [vvl] (with < v < 1) direction for the 
T, R, and Ma phases, respectively. The electrical polar- 
ization can thus be viewed as rotating in a (110) plane 
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FIG. 1: {{ux), (u y ), (u z )) Cartesian coordinates of the su- 
percell average of the local mode vectors, as a function of 
composition in Pb^Zri-zTi^Os) solid solutions. Circles cor- 
respond to disordered alloys (i.e., cti — Q2 = ctz = 0), while 
squares and losanges correspond to solid solutions exhibiting 
short-range clustering (ai = +0.2, and 02 = 03 = 0) and 
anticlustering (qi = —0.2, and 02 = 0:3 = 0), respectively. 
The filled symbols represent the z-component, while the open 
symbols are associated with the x and y components (which 
are always equal to each other). T, Ma, and R denote the 
tetragonal, monoclinic and rhombohedral phases discussed in 
the text, respectively. 



from [001] to [111] in the Ma phase as the Ti composi- 
tion decreases 0. Note that the precise compositions at 
which the T-to-MA and MA-to-R phase transitions occur 
are numerically found by identifying these compositions 

with those yielding a peak in the electromechanical re- 
let 

sponses |tj. 

Figure 1 also illustrates the effect of short-range clus- 
tering and anticlustering on the compositions delimit- 
ing the MPB, by displaying (u) as a function of x in 
Pb(Zi'i_ 2; Ti a ;)03 alloys characterized by a% = +0.2 and 
qi = —0.2, respectively. Striking features revealed by 
Fig. 1 are that SRO clustering considerably widens this 
morphotropic phase boundary, while SRO anticluster- 
ing leads to a significant compositional narrowing of the 
MPB with respect to the random case. As a matter of 
fact, samples associated with u\ = +0.2 and ot\ = —0.2 
exhibit a MPB width of 1.3% and 0.7%, respectively. 
These are quite remarkable changes with respect to the 
corresponding width of 1.0% for disordered materials, 
when recalling that the short-range orders considered 
here only occur in the first nearest neighbor shells and 
have a relatively small strength. Interestingly, Fig. 1 
shows that the T-Ma transition happens near the same 
Ti compositions of 48.7%, independently of the value and 
sign of ot\. It is thus the Ma-R transition point that is 
shifting to smaller (respectively, larger) Ti concentrations 
with increasing clustering (respectively, anticlustering). 
In other words, SRO has nearly no effect on the direc- 
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tion of the polarization in the Ti-rich area of the MPB 
of PZT, while it can considerably affect the properties of 
the Ti-poor region of this MPB. This drastic difference 
may be related to the fact that the derivative of (u) with 
respect to composition has a much smaller magnitude for 
large Ti-compositions than for small Ti-concentrations in 
the MPB of disordered PZT (see Fig. 1). Consequently 
the total electrical polarization may be less sensitive to a 
slight change of atomic configuration (i.e., a slight mod- 
ification in composition or in SRO) for larger x than for 
smaller x in the MPB of Pb(Zri_ x Ti x )0 3 . The differ- 
ent behavior of the Ti-rich and Ti-poor sides of the MPB 
may also be related to the facts that the T-M^ transition 
is second order in character while the M^-R transition is 
first order [2J| . Identifying the precise microscopic mech- 
anism responsible for the insensitivity of the Ti-rich area 
of the MPB of PZT with SRO remains for future work. 

In light of the results shown in Fig. 1, we decided to 
investigate in more detail the Ti-poor area of the MPB. 
Figure 2a shows that, as ct\ increases above zero (cluster- 
ing), the polarization of Pb(Zro.52Tio.4s)03 tends to ro- 
tate within a (110) plane and toward the [001] direction, 
since (u z ) becomes larger while (u x ) and (u y ) are reduced 
in magnitude. (We find that even for very large values of 
ot\ that are not shown here, the resulting phase remains 
monoclinic rather than becoming tetragonal.) This is 
consistent with the clustering-induced widening of the 
Ti-poor region of the MPB area displayed in Fig. 1. Fig- 
ure 2a further demonstrates that, as expected by looking 
at Fig. 1, decreasing a.\ below zero (anticlustering) leads 
to a rotation of the polarization within a (110) plane and 
toward the [111] crystallographic direction. (In fact, the 
polarization actually reaches the [111] direction charac- 
teristic of the R phase for values of a,\ below —0.2 in 
Pb(Zr .52Ti . 48)03.) Figure 2b shows the behavior of 
l/^k = 0)| 2 (sec Eq. (2)) versus a%. Figure 2c shows 
the SRO dependence of the magnitude of (u) as com- 
pared to those of (uz r ) and (uxi), the local mode av- 
erages restricted to cells occupied by Zr and Ti atoms, 
respectively. Figure 2d displays the variation of the angle 
between the pseudo-cubic [001] direction and (u), (uz r ), 
and (uxi) as a function of a\. Figure 2b reveals that 
increasing clustering leads to a larger fluctuation of the 
x- and y-components of the local polarization associated 
with the five-atoms- unit cells, as demonstrated by the de- 
crease of |/ x (k = 0)| 2 and |/ y (k = 0)| 2 as u\ increases. 
Figures 2c and 2d tell us that this large fluctuation is 
correlated with the clustering-induced enhancement of 
difference between (uz r ) and (uxi): for large values of 
ai, the local polarizations associated with cells centered 
on Zr atoms are much closer to the [111] direction — and 
much smaller in magnitude, as consistent with Refs. 
and jl^] — than the local polarizations associated with 
cells centered on Ti atoms. In other words, clustering al- 
lows the cells centered on Zr and Ti atoms to have rather 
different ferroelectric properties. 
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FIG. 2: Properties of Pb(Zro.52Tio.48)C>3 as a function of the 
Qfi SRO parameter, (a): ((u x ), (u y ), (u z )) Cartesian coordi- 
nates of the supercell average of the local mode vectors, (b): 
Fourier coefficients defined in Eq. (2) and computed at the 
r point, (c): Magnitude of the local mode vectors averaged 
over the 5-atom cells centered on Zr (circles), Ti (diamonds) 
and on all the B-atoms (solid line), (d): Same as (c) but for 
the angle between these local-modes and the [001] direction. 



Let us now discuss and understand the results depicted 
in Fig. 2. Two helpful trends can be formulated by in- 
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specting the parameters of our alloy effective Hamilto- 
nian for PZT [jj Trend I is that alloy onsite terms 
favor atomic-induced distinction between local modes, 
namely local modes associated with cells occupied by Ti 
(resp. Zr) atoms want to align along a (100) tetragonal 
(resp. (Ill) rhombohcdral) direction. Trend II is that 
the intersite interactions between a local mode cen- 
tered on the ith cell and the local modes u., associated 
with its neighbors favor similarity between and Uj . In 
other words, an homogeneous distribution of modes is en- 
ergetically favored by Trend II. When clustering of alike 
atoms occurs, Ti- and Zr-rich regions tend to develop 
dipole moments aligned along tetragonal and rhombo- 
hedral directions respectively (see Fig. 2d), in order to 
satisfy both trends I and II inside each of these regions. 
Moreover, Trend II ensures that dipole moments of these 
different regions are correlated, which results in a very 
homogeneous distribution of m, z local-mode components 
and a significantly disordered distribution of u^ x and Ui >v 
(see Fig. 2b and note that z is the preferential polariza- 
tion direction of the Ti-rich regions). Hence, clustering 
favors an overall spatially inhomogeneous monoclinic Ma 
phase in the Ti-poor area of the MPB. On the other 
hand, anticlustering leads to an incompatibility between 
trends I and II even at a very short-range scale. As a mat- 
ter of fact, local modes following their particular trend of 
being either tetragonal or rhombohcdral will find the op- 
position of their first nearest neighbors already. This re- 
sults in a more homogeneous local-mode distribution, as 
shown in Fig. 2. Very importantly, our results prove that 
this homogeneity leads to the suppression of the mono- 
clinic phase! Indeed, Figs 2 show that upon decreasing 
a% we get a distinct transition from an inhomogeneous to 
a more homogeneous distribution of local modes that is 
accompanied by a macroscopic monoclinic to rhombohe- 
dral phase transition. It is now clear why anticlustering 
results in a smaller MPB region, as shown in Fig. 1. Our 
results are consistent with Refs. and [2j|, in which 
it was pointed out that inhomogeneity is essential to the 
existence of the monoclinic phase. 

In summary, we have performed ab-initio-based simu- 
lations to study the effects of SRO — occurring over the 
first nearest neighbors' shell of the B-mixed sublattice 
— on the structural properties of PZT near its MPB. A 
competition between two different energy-induced trends 
in the Ti-poor area of the MPB leads to a considerable 
widening of the MPB area when clustering occurs, while 
yielding a significant narrowing of this phase bound- 
ary when anticlustering develops. The electromechani- 
cal properties of PZT in the Ti-poor region of its MPB 
thus drastically depend on the SRO, since they are ex- 
tremely sensitive to the direction of the polarization [j| . 
In light of these predictions, we hope that the present 
Letter will encourage the attempt of experiments aimed 
at characterizing atomic short-range order in the MPB 
of perovskite solid solutions. Such experiments are very 



frequent in metallic alloys [Tjj but (to our knowledge) 
have not yet been conducted in ferroelectric alloys. 
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